Scalar gravity: Post-Newtonian corrections via an effective field theory approach 
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The problem of motion in General Relativity has lost its academic status and become an active 
research area since the next generation of gravity wave detectors will rely upon its solution. Here 
we will show within scalar gravity, how ideas borrowed from Quantum Field Theory can be used 
to solve the problem of motion in a systematic fashion. We will concentrate in Post-Newtonian 
corrections. We will calculate the Einstein-Infeld-Hoffmann action and show how a systematic 
perturbative expansion puts strong constraints on the couplings of non-derivative interactions in 
the theory. 



Consider the apparently simple problem of the earth motion around the sun. The Newtonian solution is an 
excellent approximation but suppose we wish to be more accurate. A closer look reveals that there are many sources 
of complication. Einstein theory teaches us how to correct for relativistic effects. However, the earth is clearly not a 
point particle, and will thus deform under the influence of tidal forces. In addition the whole sun-earth system will 
radiate energy away in the form of gravitational waves. The inclusion of all these effects can make the problem of 
solving for the trajectory intractable. In the past, solving this problem was only of academic interest, but the next 
generation of gravity wave detectors will rely upon its solution [l| . The construction of accurate templates for gravity 
wave interferometers is a daunting task. After more than ten years of work the templates have been completed up to 
third Post-Newtonian (PN) order for non-spinning compact binaries 0]. However, it was not clear how to: proceed 
to higher orders in a systematic fashion, include finite size effects due to spin or spin-spin corrections. During the last 
years a new framework has emerged, coined NRGR (Non-Rclativistic General Relativity) due to its similarities with 
Effective Field Theory (EFT) ideas in particle physics, where all of the apparent obstacles of the traditional approach 
can been successfully overcome [3, |4| . NRGR naturally allows for a systematic account of the internal structure of the 
binary constitutes and permits us to calculate back reaction as well as dissipative effects Q. Moreover, new results 
for spinning compact binaries have been recently reported [fjj]. In this short contribution we will show within scalar 
gravity, how an EFT approach can be used to solve the problem of motion in a systematic fashion Q ■ In particular 
we will calculate the Einstein-Infeld-Hoffmann (EIH) action for the case of two scalar- gravitating bodies, accurate up 
to 1PN. The purpose of this contribution is pedagogical, allowing us to concentrate on the conceptual aspects. As we 
shall see a systematic perturbative expansion puts strong constraints on the couplings of non-derivative interactions 
in the theory. 



The starting point of the EFT approach consists of a theory of point particles coupled to a real scalar field <f> we 
shall call the "s-graviton" . For simplicity we will consider here a massive 3 theory in a Minkowski background, 
though we will discuss other type of models later on. The action will be given by S — S g + S pp , with 



describing the s-graviton dynamics and motion of the binary system (a — 1,2). In this equation M sets the coupling 
to matter, and A, /i, the self-interaction and s-graviton mass respectively. Also dr = y l rf- v ' dx h ' dx^ represents the 
proper time along the a-th particle and rf ,v = diag(+, — , — , — ), we work in h = c = 1 units. The choice of matter 
coupling is meant to resemble Einstein case, at least for the /loo mode, with M playing the role of the Planck Mass. 
The normalization is also chosen to mimic the graviton propagator. We could in principle add a set of higher order 
operators in the worldline action to account for finite size effects. However, tfi 3 theory is super-renormalizable and it is 
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possible to show that the n-point function is UV finite and no higher order operators are needed 1 . One other aspect 
of the super-renormalizability is the fact that a </> 3 self-interaction in four spacetime dimensions has a dimensionful 
coupling and the perturbative approach breaks down at distances of order 1/A. This is connected to IR divergences 
(in the massless limit) which appear in the perturbative expansion due to factors of X/E, with E the energy of the 
s-graviton. These IR divergences must cancel in any physical observable, such as the binding energy of the binary 
system. However, a resummation procedure is in general needed in order to achieve a finite result [8[. There are a 
few ways to overcome this. We could work in six dimensions where A is dimcnsionless, or with a IR cutoff. Instead 
we adopted a small s-graviton mass. Notice that a mass term can be produced by a tadpole mechanism, therefore a 
s-graviton mass, fi, would be naturally generated by quantum fluctuations since no symmetry prevents it. One would 
then expect \i ~ A. We will see in what follows how a well defined perturbation theory puts strong constraints in 
the self-interaction coupling of the theory. We will discuss later on under which circumstances this is a more generic 
phenomena. 



III. NRGR 



The power of the EFT formalism resides in a manifest power counting in the expansion parameter of the theory, in 
this case the relative velocity v. Here we will pinpoint the necessary steps and refer to Goldberger's contribution for 
further details Q ■ The expansion of the worldline Lagrangian leads to 
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where we have chosen x° as the worldline parameter. The propagator for the field </> appearing in L pp is still fully 
relativistic, and therefore a small velocity expansion has yet to be performed. To deal with this problem it is convenient 
to decompose the s-graviton field into potential modes (4>) with momentum scaling fc M ~ (v/r, 1/r) (notice they can 
never go on shell), and radiation modes (<E>) whose momentum scale as ~ (v/r, v/r). In the EFT spirit potential 
modes do not propagate and can be thus integrated out at each order in perturbation theory. Radiation s-gravitons 
on the other hand can appear on shell and must be kept as propagating degrees of freedom in order to reproduce the 
correct long distance physics. 



A. Power counting 



In the EFT approach one computes the effective action perturbatively, in our case in v, based on systematic power 
counting rules. In order to obtain the latter one starts with the scaling laws for the (4>, $) fields. For convenience one 
first introduces $>k, where the large momentum piece of the potential s-graviton is factored out 0, 0]- By expanding 
flU we get 

<<M*°)$ q (0)) = (2^ 3 (k + q)^") 2(k2 '^ 2) , mm) = J (2l)4 2 ( fc 2 + M 2) ^ (3) 

for the propagators. Notice that we have decided to keep the mass "non-perturbatively" to cure IR divergences, 
although we will assume fir < v in what follows, and treat it as a perturbation when allowed, in order to resemble 
the massless power counting rules and a 1/r leading order potential. If we assign the scaling x° ~ v/r wc obtain the 
following leading order power counting rules 

v 2 

4> ~ v/r, $k ~ Vvr 2 -> $ ~ M —=, (4) 

vL 

where L = mvr. The last arrow follows from the assumption that the leading order potential is given by 1/r and 
hence the virial theorem, v 2 ~ -j^r^i applies 2 . This assumption is true in the case of A = 0,/ir < v, however A / 
can easily jeopardize our power counting due to the introduction of a dimensionful coupling. In order to tackle this 



1 Notice also that using field re-definitions (~ equations of motion) we can always trade d 2 4> by a polynomial, so higher dimensional 
operator are always of the form <f> a and could be absorbed into the worldline couplings. 

2 Notice it also implies m/M ~ i/Lv 
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problem let us exercise our scaling rules and power count the first correction due to A. The diagram is shown in fig. 
2a and it scales as 



fig. 2a . 
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v 2 (XMr 2 ) times the leading order term which scales as L Q. It is easy to see higher order terms in A follow the same 
pattern. For A = g 2 M, with g a dimensionless coupling, we end up with r < -^jj for the validity of the perturbative 
approximation and power counting. In order to make sense of the perturbative approach we had to cure this IR 
singularities before expecting any power counting to work, and that is what the s-graviton mass is doing. If we 
demand our leading order potential to match the Newtonian case we will set M = mpi and therefore the perturbative 
expansion is valid for r < Ipi/g, with Ipi = 10 -33 cm, the Planck length. To avoid entering the quantum realm we 
will have to fine tune g to an extremely small number of the order of 10 -40 for typical binary systems in the solar mass 
range in the inspiral regime. This obviously defies naturalness arguments and puts a flag on the phenomenological 
viability of such theory since it implies a ridiculously small self-coupling, A ~ 10~ so mpi ~ 10~ 60 GeV\ Notice that the 
problem does not lie in the self-coupling itself but in the strength of the worldline coupling which determines the leading 
order scaling laws. In Einstein theory this is taken care of by the three graviton coupling, g^ ~ k 2 /M — ► g^Mr 2 ~ 1. 
The condition fj, < v/r also implies a stunningly tiny s-graviton mass of the order of 10~ 30 ev. These are consistent, 
and somehow equivalent, to solar system constraints Q, whereas by naturalness arguments fi ~ A would produce an 
even smaller value. We will hereon assume \ir ~ v, Xr -C 1, and proceed with this theory as a playground. 



IV. EINSTEIN-INFELD-HOFFMANN 



Let us concentrate now in the calculation of the 1PN correction to the gravitational potential. The leading order 
one s-graviton exchange can be easily seen to reproduce Newtonian gravity Q. We also need to take into account 
diagrams with one single s-graviton exchange which are down by v 2 shown in figures la and lb plus the non linear 
terms depicted in figures 2a and 2b. We will treat the s-graviton mass as perturbation in the one s-graviton exchange 
(pr ~ v) and that is shown in diagram lc. The computation proceeds systematically by using the Feynman rules of 
the EFT order by order. We will concentrate in detail in fig. 2a, we will display the full result later on. For the three 
s-graviton diagram we will have 

fig. 2a = 1 / dMMMn^iM^)^')))- (6) 
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Our task now is to compute the three-point function. For a <fi 3 theory one obtains 

d 3 k r 



(TWxiWxt^M)) = 3l(-iX)S(t 1 -t 2 )S(t 1 -t 2l ) fU^e-^^pirfS* II 2(k ? + (7) 

The next step would be to plug this expression back into © , get a finite result which we will have to further expand 
in powers of \ir ~ v and keep the leading order piece, already at 1PN for XMr 2 ~ 1. In the EFT spirit a better way 
to proceed is to treat /x as a perturbation in the same way time derivatives are treated, by expanding the propagators 
in powers of /x/|k|. For the one s-graviton exchange this represents no harm. In general one faces the problem that 
IR divergences will only cancel out after all the terms are included. If we are willing to accept that is the case one can 
calculate the 1PN correction by taking the massless limit of and keep the (non-constant) finite piece. Therefore, 
introducing d = 3 + e and taking the limit e — > one gets 

fie 2a - i} 3mim " fdt d3k2d3kl 1 c-ikx-fxx-x,) . . iX ^G N m im 2 f A, 1 kl . (xi _ X2) 

ng. Za - IA ^3 J ar (27r) 3( 27r )3 k 2 k 2 (ki+k2) 2 e ~ %A lgM J M ( 27r )3 ( k 2 )3/2 

^^r(e/2) [dt( \ Xl ~ X2 \ \ — > -iSXMG 2 N m\mi [ dt log(dxi - x 2 |) + constant. (8) 
■KM J \ 4 J J 



ax 



3 Recall that loop effects in NRGR for gravitons are suppressed by 1/L and can be thus ignored in the classical scenario Q . 
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FIG. 1: One s-graviton exchange contribution at 1PN. The (^) represents a correction to the propagator, and x a mass insertion. 
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FIG. 2: Non linear contributions at 1PN. 

with Gn = M' 1 > ana - ^he "constant" piece also contains the \ IR pole 4 . 

Our final task consists in collecting the other few pieces. We refer to [1, 0| for details since the calculations are 
almost identical. Let us compute the result for the new term in fig. lc due to the s-graviton mass insertion, 

fig. lc = -i^^ J dhdhSih - t 2 ) J j£ e -*-(~-*0 = 1 j dt G N m im2 ^\ Xl - x 2 |, (9) 



which is nothing but the 0(v 2 ) piece in the expansion of the Yukawa potential, —^-p— ~ /i — -(1 + ^ — l~ •••)• 
Putting everything together, including mirror images, we finally obtain 



Leih = lj2 maV * 



2|xi -x 2 | 



2 , / \ ( V l ' X1 2 )(V 2 • Xi 2 ) 
V 2 + (Vl • V 2 ) —2 



Xl - x 2 



G 2 N mim 2 (mi + m 2 ) 
2|x 1 -x 2 | 2 



+ 2 ^ m i m 2M 2 |xi - x 2 | - 3A G^r M m 2 mi(mi + to 2 ) log(/i|xi - x 2 |) (10) 

where we have also included the relativistic corrections to the kinetic energy of the point particles. The logarithmic 
potential introduces a very interesting feature, namely a - force into the equations of motion and therefore v 2 ~ 
a/r + b + which implies a dark matter type of effect for the galaxy rotation curves. This is however by no means 
a serious candidate and we mention this only as a curiosity. 



V. DISCUSSION - CONCLUSIONS 



The EFT approach is a powerful tool within the PN framework. Using no more than dimensional analysis many 
conclusions can be already drawn before dwelling into the details of the calculations. We applied the techniques in 
the case of a massive 3 theory as a playground but the ideas can be easily extended to more complicated scenarios. 
From the NRGR power counting rules we learned that in order to produce a well defined perturbative expansion, A 
had to be fine tuned to a ridiculously small (compared with mpi) scale. One might however ask whether this is a 
feature of a <fi 3 theory or will this be faced in other scenarios. Let us consider a more general case, 

J d A x ( 7 (0)c^cKV + B(4>)) (11) 

In what follows we will consider two distinct case. 



4 In the massive case this is translated into a log/i factor. The full result is a Bessel function, Ko(fir), whose leading order piece in [ir 
reproduces the logarithmic potential. 
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A. B{(j}) = 

If wc expand ^{(j)) ~ 1 + (f>/M + we will get a kinetic piece plus a potential V(<p) with terms like (4f \ <t>d 2 (j) 

n > 1. This theory is not renormalizable, and it is easy to show it resembles Einstein case. We can also show that 
the pcrturbative approach is under control by power counting the contribution from a generic term in V(cf>). For an 
(n + 2) s-graviton diagram we will get, 

/ 2 \ ™+ 2 

VZ" + -L= ) M 2 r 2 /v ~ i« 2 ". (12) 



For instance the first term in the expansion, </> 2 <9 2 0, resembles the three graviton coupling in Einstein theory (up 
to tensor structure). Had we chosen this interaction we would have ended up with a similar 1PN correction as in the 
original Einstein-Infeld-Hoffmann action [H, 0] • 



B. B(4>)^0 

This case is substantially different. Let us study a generic term, gM ((f)/M) n , with g a dimensionless coupling and 
n > 4. The n s-graviton diagram will scale as 

jMV 4 !) 2 "" 1 ~ gL 2 v 2n ~ 7 ~ L g(m/M) 2 v 2n - 8 . (13) 

To have a controlled perturbative expansion wc would have to impose gv 2n ~ 8 (m / M) 2 < 1. For the marginal case 

n = 4, setting M = rnp; one needs g < (^^p-) ~ 10 -70 ! for solar mass binary constitutes. Wc can improve this 
number by considering higher dimensional terms, namely larger n, but the enhancement is really minute. Notice that 
this problem arises at the classical level since the coupling to elementary particles is already too small to represent 
any trouble. Is only in the superposition of terms, which build up the massive lump of the star, that the perturbative 
expansion breaks down. From this analysis we conclude that in pure scalar gravity non-derivative self-interactions are 
extremely constrained. 

One could then wonder about more general models including scalar fields, like tensor-scalar gravity [1(3]. In the latter 
in addition to the graviton field a scalar interaction is added with an action similar to pip in a curved spacctimc 
background. Within this type of scenarios the problems we encountered here can be cured by modifying the power 
counting. For instance, a large mass can be added to the scalar field (larger than the inverse of the solar system 
distance), which will render the field a negligible short range interaction. Another possibility would be to keep it 
nearly massless but weaken the coupling to matter to a much feeble strength M ^> mpi. In this case the 3-scalar 
diagram (fig 2a) will now scale as (jj) 3 ■ For A ~ gM one needs g{m/M) 2 < v 2 in order to be treat as a 

perturbation. By naturalness argument one would expect g ~ 1, and we will then have a very tiny coupling to 
elementary particles. For instance, the coupling to a proton will be of the order of m prot on/M ~ 10 -60 ! For a 4 

theory the 4-scalar diagram would now scale as ^(m/M) 4 . Compared to the leading Newtonian potential we get a 

suppression of order \(m/M)' 3 (Mrv 2 )~ 1 , which can be seen to be effectively small for A ~ 1, M ~ m. Again the 
coupling to elementary particles is highly suppressed. Both solutions will keep the theory consistent with experimental 
data, both rely however in the introduction of a high mass scale into the theory, much higher than the Planck scale 
or the scale of particle physics. Perhaps this is an indication that non-derivative self-interactions are not present in 
nature 5 
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5 The possibility that the dilaton, the scalar Goldstone mode of the SSB of conformal invariance, could provide an analog of Einstein 
theory was recently raised in fllll . Contrary to the traditional lore shift invariance does not forbid a non-derivative interaction as in the 
case of internal symmetries. A similar fine tunning, which was argued in to be of the same nature as the cosmological constant, was 
invoked to avoid unstable configurations. 
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